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Abstract 



Assembly and disassembly dynamics of microtubules (MTs) is tightly controlled by 
MT associated proteins. Here, we investigate how plus-end-directed depolymerases 
of the kinesin-8 family regulate MT depolymerization dynamics. Employing an 
individual-based model, we reproduce experimental findings. Moreover, crowding 
is identified as the key regulatory mechanism of depolymerization dynamics. Our 
analysis gives two qualitatively distinct regimes. For motor densities above a par- 
ticular threshold, a macroscopic traffic jam emerges at the plus-end and the MT 
dynamics become independent of the motor concentration. Below this threshold, 
microscopic traffic jams at the tip arise which cancel out the effect of the de- 
polymerization kinetics such that the depolymerization speed is solely determined 
by the motor density. Because this density changes over the MT length, length- 
dependent regulation is possible. Remarkably, motor cooperativity does not affect 
the depolymerization speed but only the end-residence time of depolymerases. 

Key words: kinesin-8; length-regulation; microtubule dynamics; traffic jam; 
particle exclusion; driven transport 



INTRODUCTION 



Microtubules (MTs) are cytoskeletal filaments that serve a central role in intra- 
cellular organization (JTJ [2j) and several cellular processes including mitosis ([31 2]), 
cytokinesis (j5j) and intracellular transport (6). They can cope with this multitude 
of diverse tasks because they are highly dynamic structures which continually as- 
semble and disassemble through the addition and removal of tubulin heterodimers 
at their ends. GTP-hydrolysis is the energy source which drives switching between 
persistent states of growth and shrinkage, a stochastic process termed dynamic in- 
stability (IT HTOj) . Each cellular process employs a specific set of MT-associated 
proteins (M APs) to tightly regulate the rates of growth and shrinkage as well as 
the rate of transition between these states (fTTHT3"j) . 

Depolymerases from the kinesin-8 and kinesin-13 protein families (e.g., Kip3p 
and MCAK, respectively) are important regulators of MT dynamics. They are 
thought to promote switching of MTs from growth to shrinkage (catastrophes) 
(fT2"j) . Whereas MCAK lacks directed motility and diffuses along MTs ( I14p . Kip3p 
is a highly processive plus-end-directed motor (|l5| H6 p . Proteins from the kinesin-8 
family are important for regulating MT dynamics in diverse organisms. Kifl8A 
is a key component in chromosome positioning in mammalian cells (1TT1 - [T9]) where 
it regulates plus-end dynamics. Its orthologs, the plus-end directed motors Kip3p 
in budding yeast (1161) and Klp5/6 in fission yeast ( l2Ul - [2"2"]) . show depolymerizing 
activity. A notable feature shared by these MT plus-end depolymerases is that they 
depolymerize longer MTs more rapidly than they do shorter ones ( 1T51 ITT1 |2T| 1231 . 
A similar length-dependent regulation of MT assembly by kinesin-5 motors was 
observed in in vivo studies of chromosome congression in budding yeast ( 12^|) . The 
key experimental observations from in vitro studies of Kip3p (23) are that 1), the 
end-residence time of Kip3p at the tip depends on the bulk concentration of Kip3p 
and correlates inversely with the macroscopic depolymerization speed; and 2), the 
macroscopic depolymerization rate is directly proportional to the flux of Kip3p 
towards the MT plus-end. 

It is thought that length-dependent depolymerization kinetics serves several 
purposes §Z§. For example, positioning of the nucleus at the cell center during 
interphase is achieved by growing MTs that push against the cell poles while re- 
maining attached to the nucleus. A higher rate of catastrophes for longer MTs 
implies that shorter MTs have an increased contact time with the cell poles. Com- 
puter simulations show that this leads to a higher efficiency of nuclear positioning 
during interphase (I25p . 

There is convincing experimental evidence that molecular traffic along MTs 
strongly affects MT depolymerization dynamics. However, in vitro experiments 
can not yet fully explore the underlying traffic dynamics. Theoretical investigations 
employing individual-based models can be instrumental in furthering a mechanistic 
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understanding of this process. Fortunately, these models can be constructed on the 
basis of substantial quantitative data available from in vitro experiments ([T5] [25]) 
characterizing the binding kinetics and the motor activity of plus-end-directed 
motors. Therefore we sought to identify the molecular mechanisms underlying 
the observed correlation between depolymerization dynamics and molecular traffic 
along MTs. 

In this study, we constructed an individual-based model for the coupled dynam- 
ics of MT depolymerization and molecular traffic of plus-end-directed motors. This 
model quantitatively reproduces previous experimental results flT5 | 123]) . Moreover, 
we make precise quantitative predictions for the density profiles of molecular mo- 
tors on the MT and demonstrate that molecular crowding and ensuing traffic jams 
regulate the depolymerization dynamics. We find two qualitatively distinct regimes 
of depolymerization dynamics: At low bulk concentrations of depolymerases, the 
depolymerization speed of MTs is density-limited and is a function of the bulk 
concentration and average motor speed alone. There is a sharp threshold in bulk 
depolymerase concentration above which macroscopic traffic jams emerge and the 
depolymerization speed is simply given by the microscopic depolymerization rate. 
Of note, none of these features are affected by the degree of cooperativity in the 
depolymerization kinetics. In contrast, the end-residence time of a depolymerase 
(i.e., the typical time it spends at the plus-end) is strongly correlated with coop- 
erativity. We outline how these predictions from our theoretical analysis can be 
tested experimentally. 

RESULTS 

Model definition 

We use an individual-based model, as illustrated in Fig. [TJ to describe the dy- 
namics of plus-end-directed depolymerases. Motor proteins, present at a constant 
bulk concentration c, are assumed to randomly bind to and unbind from the MT 
lattice with rates u a and Ud, respectively. Bound motors are described as Pois- 
son steppers (A more detailed biochemical model for motors on microtubules has 
to await further experimental analysis. One of the different possible schemes has 
recently been studied by Klumpp et al. (26).) that processively walk along individ- 
ual protofilaments towards the plus-end at an average speed v (1271) . These motors 
hinder each other sterically because individual binding sites % = 1, . . . , L on each 
protofilament can either be empty (rij = 0) or occupied by a single motor (rij = 1). 
Because switching between protofilaments is rare (127]) . transport along each of 
the protofilaments can be taken as independent, and the model becomes effec- 
tively one- dimensional (12"%]) (Fig. ITlB) . Models of this type were recently discussed 



3 




FIGURE 1: Illustration of MT and motor dynamics. Molecular motors present at 
concentration c randomly attach to unoccupied tubulin dimers along the MT lat- 
tice with rate u a . While bound they processively move toward the plus-end at rate 
v, and unbind with rate u^,. Because motors do not switch lanes (protofilaments) , 
the MT lattice (A) becomes effectively one-dimensional (B). Each lattice site 
(with i — 1, . . . , L numbering the sites) may be empty (n^ = 0) or occupied by a 
single motor (n, = 1). At the plus-end the motors act as depolymerases (indicated 
by scissors) either alone with rate So or cooperatively with rate Si. 

as minimal models for intracellular transport In its given formulation, 

where the cytosol is considered as a homogeneous and constant reservoir of motors, 
it is equivalent to a driven lattice gas model known as the totally asymmetric sim- 
ple exclusion process with Langmuir kinetics (TASEP/LK) ( 129]) . A central finding 
of this model is that the interplay between on-off (Langmuir) kinetics and directed 
transport along protofilaments can result in "traffic jams" in which the density pro- 
file of motors along a protofilament shows a sharp increase from a low-density to a 
crowded high-density regime (1291 13TI) . Such and other crowding effects (1331 l34"l) are 
important for a molecular understanding of MT dynamics. Previous theoretical 
studies on this topic largely disregarded crowding effects or considered parame- 
ter regimes in which they are unimportant (p5oT[571) . Depolymerization, including 
crowding effects, has also been investigated for diffusive depolymerases such as 
MCAK (EE]). 
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At the plus-end of the system we consider depolymerization dynamics arising 
due to the interaction of molecular motors with the MT tip. Motivated by recent 
experiments (123]) . we assume non-processive depolymerization, i.e, a molecular 
motor dissociates from the lattice after triggering depolymerization. Because the 
molecular mechanisms are not yet fully resolved, we study two scenarios of de- 
polymerization (see Fig. In the noncooperative scenario, the dissociation rate 
depends only on whether the last site is occupied by a motor. If the last site is oc- 
cupied, m — 1, the MT depolymerizes at rate do- However, recent single molecule 
studies indicate that Kip3p may act cooperatively fl23j) . which we consider as our 
second scenario. After arriving at the plus-end, the motor is observed to pause 
and depolymerize a tubulin dimer only after a second Kip3p has arrived behind 
it. In this scenario, a tubulin dimer is depolymerized with rate 5\ if both the last 
and the second-to-last sites are occupied, rii-i — Ul = 1. Therefore, the total 
depolymerization rate can be written as: 

A = 5 n L + b~in L _in L . (1) 

For stabilized MTs, the spontaneous depolymerization rate is small (123]) and thus is 
not considered here. The relative magnitude of the noncooperative rate 5q and the 
cooperative rate S\ determines the degree of cooperativity of the depolymerization 
kinetics. In an average over many realizations of the stochastic process (ensemble 
average), the depolymerization speed Vd epo i depends on the occupation of the last 
two binding sites by depolymerases (Fig. [II 



Vdepoi = (o~oP+ + SiK + ) a , (2) 

where a is the lattice spacing. Here p + := (til) is the probability that the last site 
is occupied (i.e., the expected motor density at the plus-end), and k + := (t^l-i^l) 
denotes the probability that both the last and second-to-last sites are occupied. 
We analyzed this model via stochastic simulations and analytic calculations (for 
further details, see the Supporting Material). 



Validation of the model and its parameters 

The model parameters far as they are available, fixed by experimental data. 

The motor speed, v , the motor run length, £, and motor association rate, u a , were 
measured previously (1231) : 

v = 3.2 yummin -1 , 
u a = 24nM -1 min -1 / um~ 1 , 
£ m 11 jjm . 
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Using an MT lattice spacing of a = 8.4 nm, we derive the corresponding parameters 
in our model as follows: The motor speed v corresponds to 6.35 lattice sites per 
second, i.e., a hopping rate of v = v/a = 6.35 s -1 . The inverse hopping rate 
r := v~ x = 0.16 s and the size a of a tubulin dimer serve as our basic timescale 
and length scale, respectively. Then, the measured association rate corresponds 
to a rate u a ~ 5.3 x 10~ 4 nM _1 site _1 r _1 . The dissociation rate, Ud = v/£, is 
derived as the ratio of the mean motor speed, v, and the mean motor run-length, 
£. The latter equals 1310 lattice sites. Thus, the dissociation rate is expressed 
as Ud ~ 7.6 x 10~ 4 site -1 r _1 . In contrast to the transport behavior on the MT, 
the parameters concerning the depolymerization rates, <5q/i, cannot be directly 
extracted from experiments. However, there is evidence for a depolymerization 
rate as high as the motor speed, v [231 . As a starting point for the following 
discussion we tentatively take So = v. 

Using the above set of parameters we now phenomenologically compare the 
results from numerical simulations of our model with observations from experi- 
ments. Specifically, we consider kymographs of the MT, which show how the MT 
length and the motor density on the MT evolve over time. For the simulation data 
shown in Fig. [2] we consider an MT consisting of 14 independent protofilaments 
and investigate the dynamics for the noncooperative scenario and a range of motor 
concentrations, c = 1.2, 1.8, 2.6 nM, cf. Fig. [2]A-C. Surprisingly, as shown later, 
neither the cooperativity of the motors nor a decrease in the depolymerization 
rates led to different shapes of kymographs (see also Fig. SI). We find an initial 
time period in which, starting from an empty MT lattice, the motors first fill up 
the lattice ( )39| l4"0l) . This is followed by a time window in which the motor den- 
sity exhibits a quasi-stationary profile, i.e., the density at a certain distance from 
the minus-end does not change except for boundary effects induced by the plus- 
end. The corresponding density profiles are illustrated in Fig. [^i? and discussed 
in more detail in the following section. In this quasi-stationary regime, the de- 
polymerization dynamics shows qualitatively different behavior depending on the 
concentration of free motor molecules: At low concentration, c < 1.4 nM, and thus 
low density of motors on the MT, depolymerization slows down gradually in the 
course of time (Fig. |2]/1). When the motor concentration increases to larger values, 
c > 1.4 nM, an intermediate regime emerges in which the depolymerization speed 
stays roughly constant (Fig. [2]B and O). Remarkably, we find that during this 
regime, the depolymerization speed is directly proportional to the motor density, 
Vdepoi(£) = P-{L)v (Fig. [2]D). At a third stage in the depolymerization process, 
there is a rather abrupt change in the depolymerization speed right where the 
density profile shows a steep drop (Fig. [2JC-.E). After we have elaborated more on 
the theoretical model, we discuss why there is such a tight correlation between the 
depolymerization dynamics and the density profile. 



6 




Time t [min] MT length L [/mi] Kip3p concentration c [nM] 



FIGURE 2: Validation of the theoretical model. (A)-(C) Time-space plots of 
stochastic simulations for a range of motor concentrations and depolymerization 
rate 5q = 6.35 sites s^ 1 . The density of molecular motors is shown as the bright 
area (green), and the MT is shown as the dim area (red; for details, see Supporting 
Material). For low concentrations, c < 1.4 nM, depolymerization slows down grad- 
ually (23). At higher concentrations, c > 1.4 nM, there is a rather abrupt change in 
MT shortening. This change is correlated with a steep decrease in the motor den- 
sity (DW), indicated as dotted lines. (D) The depolymerization speed, Vdepoi ; as a 
function of the length of the MT L(t), extracted from the simulation data shown 
in the kymograph (gray). The position of the DW (dotted), and the predicted 
depolymerization speed, Vdepoi — v p(L) (see also Eq. 10), using the linear approxi- 
mation for the motor density profile (black) and the density profile extracted from 
stochastic simulations (green), coincide very well with the observed depolymer- 
ization speed; v = 6.35 sites s _1 is the walking speed of the motors. (E) Density 
profiles at the minus-end from stochastic simulations (lines with symbols), exact 
solutions (solid) and linearized theory (dotted) are shown. (F) As a function of 
the motor concentration, c, and the distance from the minus-end, there are dis- 
tinct types of density profiles. At motor concentration lower than c = 1.4 nM 
(thin black), the density of motors along the MT is low and the profile is smooth. 
The Langmuir density is reached continuously after a certain MT length (dashed, 
numerical). At high concentrations, c > 1.4 nM there are two regions along the 
MT separated by an intervening DW (black, exact; see Supporting Material): an 
approximately linear antenna profile and a flat profile (Langmuir density). Linear 
approximations for the continuous and the discontinuous transitions (Eq. [4]) are 
shown as well (gray). Thin lines refer to the density profiles shown in E. 
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All of these qualitative features of MT dynamics are identical to those found ex- 
perimentally ([TBI I2"3"j) . and suggest that the density profile and, in particular, traffic 
jams formed on the MT lattice are the main determinants of the depolymerization 
dynamics. Moreover, the time scales of the dynamics agree quantitatively well 
with experimental results for the same motor concentrations ( JT5| [23]) . This vali- 
dates our theoretical model because up to the depolymerization rate 5, all of the 
model parameters were derived from experimental data (I23p . 



Density profiles at the minus-end (bulk density) 

The above observations strongly point toward a tight correlation between the de- 
polymerization speed and the motor density profile at the minus end, p-(x), which 
we henceforth call the bulk (motor) density. The quasi-stationary bulk density pro- 
files shown in Fig. [2]E were obtained by assuming very long lattices; effects caused 
by the plus-end are not visible in the vicinity of the minus-end. A more detailed 
discussion of these simulations can be found in the Supporting Material. Since this 
bulk density will play an important role in the following analysis, we summarize 
its features as obtained from analytical calculations detailed in the Supporting 
Material. 

At the minus-end, the density profiles show an initial linear increase. This is 
an "antenna effect" (15) as illustrated in Fig. [3]/L Motors that attach in proximity 
of the MT minus-end immediately move toward the plus-end, thereby generating 
an approximately linearly increasing accumulation of motors. The slope is given 
by K/£, where K = cu a /ud denotes the binding constant. At sufficiently large 
distances from the minus-end, the density profile becomes flat and dominated by 
Langmuir kinetics with the ensuing Langmuir density: 

K cu a 

PLa = : , r? = , • ( 3 ) 

1 + K CU a + U d 

The full density profile is obtained by concatenating the antenna profile and the 
flat Langmuir profile such that the motor current is continuous along the MT. 
We find two qualitatively distinct scenarios (Fig. j^Jfi 1 ). For low concentrations of 
molecular motors, c, the antenna profile matches the asymptotic Langmuir density 
continuously, resulting in a wedge-like profile. In contrast, above a certain thresh- 
old value for the concentration, determined by the binding constant K~ = 1, the 
two profiles can no longer be matched continuously and the density profile displays 
a sharp discontinuity, also termed a "domain wall" (DW) (129|) . In other words, if 
the Langmuir density rises above a critical value of p£ a = 0.5, a crowding-induced 
traffic jam will result (jUJ) (Fig. [3]A). The density profiles obtained from the an- 
alytic calculations and the stochastic simulations agree nicely, as illustrated in 
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FIGURE 3: Illustration of the antenna and crowding regimes, and of cooperativity. 

(A) Starting from an empty MT, motors start to accumulate on the MT lattice 
by attachment and subsequent transport to the plus-end. The combined effect of 
Langmuir kinetics and steric exclusion between the motors leads to two sharply 
separated regimes. Starting from the minus-end, the motor density increases lin- 
early (antenna profile). At a certain critical length i~, a macroscopic traffic jam 
arises because particles hinder each other and crowding dominates the MT density. 

(B) Illustration of non-cooperative (B, nc) and fully cooperative (C, fc) depoly- 
merization kinetics. With regard to the the depolymerization speed, both models 
are effectively equal (see main text). 



Fig. (2j£'. In particular, the theoretical analysis gives an explicit expression for the 
width of the antenna-like profile: 

i+k iovK<l, 
koTk) forK>l. 

This result reduces to the average run length of molecular motors, £ = v/ud, 
in the limit of very low binding constant, K -C 1, where crowding effects can be 
neglected (1371) . However, with increasing K the regime with an antenna-like profile 
becomes significantly shorter than i (Fig. [2]F) . 
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Depolymerization dynamics is independent of cooperativity 



We now address how the cooperativity of the depolymerization kinetics affects the 
macroscopic depolymerization speed. There are two limiting cases: noncooperative 
depolymerization (nc) with (So, Si) = (6,0), and fully cooperative depolymeriza- 
tion (fc) with (So, Si) = (0, S) (for an illustration, see Fig. |3j .Band C). Remarkably, 
we find from our stochastic simulations, shown in Fig. |4| that there is no difference 
in depolymerization speed for these two limiting cases. Even when the depolymer- 
ization dynamics contains cooperative as well as noncooperative terms, we do not 
find any significant differences in the depolymerization speed (Fig. |4]Z?). 

This observation from our stochastic simulations can be explained by the follow- 
ing molecular mechanism: Consider a model with fully cooperative depolymeriza- 
tion kinetics. Then, after the first motor has arrived at the plus-end, the terminal 
site of the MT will remain occupied from that time on. Depolymerization only 
occurs if another motor arrives at the second-to-last site. In other words, while the 
last site remains occupied, the second-to-last site triggers the depolymerization. 
Hence, as far as the depolymerization speed is concerned, the fully cooperative 
model is identical to a noncooperative model with the same molecular rate S. In 
the noncooperative model the terminal tubulin dimer is removed at rate 5 once a 
molecular motor has arrived at the last site (see Fig. [3]B). In the fully cooperative 
model, the terminal tubulin dimer is removed once a molecular motor has arrived 
at the second-to-last site next to a permanently occupied last site (Fig. [3]C). 



Depolymerization dynamics is strongly affected by crowding 

To gain further insights in the correlation between the depolymerization speed and 
the density of motors on the MT, we performed stochastic simulations focusing 
on the MT plus-end by regarding the dynamics in a co-moving frame. Instead 
of simulating the full-length MT with an antenna profile and a subsequent flat 
Langmuir density, we considered a reduced model in which the density at the left 
end is set equal to the Langmuir density pLa- For long MTs, the Langmuir density 
is always reached, so that the reduced system is fully equivalent to the original 
model. Our simulations show two clearly distinct regimes of depolymerization 
dynamics (Fig. [4]): For small microscopic depolymerization rates, St < p^, the 
depolymerization speed is rate-limited: Vdepoi — a S- In contrast, for rates St > 
PLa, the depolymerization speed is density-limited, and the Langmuir density is 
the limiting factor: Vdepoi — Pl&v. The boundary between the two regimes is 
remarkably sharp and given by 

PL = St . (5) 
This implies that the depolymerization speed can switch between being density- 
limited and rate-limited by changing the concentration c or the values of the bio- 
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chemical rates of depolymerases binding to and unbinding from the MT lattice. 
Overall, the depolymerization speed obeys a scaling law 

I f° r $ T < PLa 

Vdepol = PhaV V(dT/> La ) = < r r > ( 6 ) 

[p Lsi v for dr > p La 

where V(x) is a universal scaling function with the simple form V(x) = x for x < 1 
and V(x) = 1 for x > 1. Experimentally, this implies that one should find data 
collapse upon using such a scaling plot (Fig. Ill 



To gain a molecular understanding of these remarkable features of the depoly- 
merization speed, one needs to have a closer look at the density profile of the 
molecular motors at the MT tip. If the depolymerization rate is small, 8 < u, 
motors leave the tip more slowly than they arrive. Therefore, the MT tip acts as a 
bottleneck for molecular transport that disturbs the density profiles either locally 
or macroscopically. A weak bottleneck induces a local perturbation ("spike") (133|) . 
These spikes are sharp changes of the density profile with a typical extension that 
scales with the size of a heterodimer. However, if the strength of a bottleneck 
exceeds a threshold value, the spike extends to a macroscopic perturbation ("traffic 
jam") (j33j) . Fig. |5]A illustrates how, for a given Langmuir density, pL a = 2/3, the 
effect on the density profile changes from a spike (blue) to an extended traffic jam 
(red and green) when the depolymerization rate is 5. 

Let us now analyze the conditions and consequences of such bottlenecks in 
more detail. Suppose we are in a parameter regime where the plus-end disturbs 
the density profile only locally, i.e., on the scale of a heterodimer. Then, we may 
take the bulk density to be equal to the Langmuir density, pL a , up to the last site 
(the plus-end) where it jumps to some higher or lower value p + . The particle loss 
current at the plus-end due to MT depolymerization is then given by 

^depol = (1 - PLa)p+<$ • (7) 

The factor 1 — pL a arises because the particle number decreases only if a particle 
depolymerizes the MT and the second-to-last site, L — l, is unoccupied. Otherwise, 
depolymerization dynamics and the associated frame shift of the MT lattice do 
not change the occupation of the last site. This particle loss has to be balanced 
by the incoming particle flux, 

^La = PLa(l - PLa)^- (8) 

Equating these particle fluxes (Eqs. [7] and [8]) implies the following condition for 
the motor density at the plus-end: 

(plu/St forp La <(5r 

P+ = \ -, , . r ' ( 9 ) 



1 for p La > 5t ' 
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FIGURE 4: Scaling plot for the depolymerization speed Vdepol- (A) Upon rescaling, 
both the macroscopic depolymerization speed, Vdepol , and the microscopic depoly- 
merization rate, 5, with the Langmuir density, p^ a , all data collapse onto one uni- 
versal scaling function V (solid gray). A sharp transition at 5t = p^ a distinguishes 
the rate-limited regime from the density-limited regime. (B) Comparison of co- 
operative and noncooperative depolymerization, with the macroscopic depolymer- 
ization speed Vdepoi as a function of Langmuir density pLa- For 5 := 5o + 5± = 0.7 v 
different degrees of cooperativity are displayed as indicated in the graph. 
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FIGURE 5: Density profiles at the plus-end, corresponding phase diagram, and 
depolymerization scenarios. (A) Density profiles at the MT plus-end in the co- 
moving frame for c = 2.9nM, and 5 = 0.1,0.3 (left), 0.35,0.5 (middle) and 0.8v 
(right). The simulation results and analytical solutions (black; see Supporting Ma- 
terial) agree nicely. (B) Depending on the value of 5 and the density of motors, 
PLa, there are three different classes of density profiles at the plus-end: wedge-like 
(diamonds) , traffic jams with a DW (square), and spikes (circles). The transition 
between profiles with an extended traffic jam and a localized spike (solid line) also 
marks a qualitative change in the depolymerization speed. Whereas the depoly- 
merization speed is density-limited in the spike regime, it is rate-limited in the 
DW and wedge regime. Symbols correspond to parameters as displayed in panel 
A. (C) Depending on the value of S and the density of motors, p^ a , there are three 
different regimes of depolymerization dynamics. In regime a depolymerization is 
density-limited for arbitrary MT length. In contrast, depolymerization is rate- 
limited for long MTs and density limited for short MTs in the regimes (3 and 7. 
For details, see the main text. 
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where the fact that the motor density is bounded p + < 1 is already accounted 
for. The particle density on the last site, in turn, determines the depolymerization 
speed. For pL a < St, one obtains according to Eq. [2] and Eq. [9] 

Vkepoi = p+Sa = p La v . (10) 

Remarkably, here the effect of the depolymerization kinetics (5) cancels out such 
that the macroscopic depolymerization speed is independent of the molecular de- 
tails of depolymerization kinetics and solely determined by the Langmuir density, 
i.e., the motor density in the bulk, p_(x), and not at the tip of the MT. This 
result crucially depends on the presence of a microscopic spike. It explains the 
hitherto puzzling experimental result that the depolymerization speed is directly 
proportional to the bulk motor current along the MT ( 12 3 p (Fig. S2). 

Because the density is bounded, p + < 1, density profiles with a spike are 
only possible if the densities are not too large, pL a < St. This is the case for 
the blue curve in Fig. [5]/L For densities exceeding the critical density, p£ a = St, 
the bottleneck-induced perturbation in the density profile can no longer remain 
a local spike, but has to become macroscopic in extent (33) (see green and red 
curves in Fig. [5]A and Supporting Material). One finds that over an extended 
region, the binding sites at the plus-end then remain permanently occupied such 
that p + = 1. This immediately implies that the depolymerization speed becomes 
density-independent and proportional to the microscopic depolymerization rate: 

Vde P oi = a<5- (11) 

There is a tight correlation between the shape of the density profiles and the 
macroscopic depolymerization speed. The analytic results explain the molecular 
mechanism behind the numerically observed scaling law (Eq. [6]), with a sharp 
transition from density-regulated to a rate-limited depolymerization dynamics at 
a critical value of p£ a = St (cf. the classification of density profiles and depoly- 
merization regimes shown in Fig. |5]S). 

Actually, the above calculations can be generalized to the regime in which the 
motor density exhibits an antenna-like linear profile, i.e., for MT length shorter 
than t~ . As detailed in the Supporting Material, we find that the depolymerization 
speed is rate-limited, Vd epo i = a,5, if MTs are shorter than t~ but still longer than 
a second threshold length: 

£ d := 5a/cuj a = £6t/K. (12) 

In contrast, for £ d > £~ , the depolymerization speed in the antenna regime is always 
length-dependent and strictly follows the shape of the antenna profile, P-(x): 

Vkepoi = P-(L)v. (13) 
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Using Eq. |4j the condition l& > t~ on the threshold lengths is equivalent to 
5t > for K < 1 and to 5t > 1 — p La for K > 1. 

Combining all of the above results, we find three mechanisms governing de- 
polymerization dynamics, as illustrated in Fig. [5]C: 

(a) For 5t > pLa, the depolymerization speed is always density-regulated and 
given by Vd epo i{L) = P-(L)v, where L is the time-dependent length of the 
MT. In this parameter regime, the depolymerization speed is a direct map 
of the bulk motor density profile on the MT, p~{x), a feature that can be 
exploited experimentally to measure the profile. 

For pLa > St > 1 — pLa, the depolymerization speed is rate-limited for MTs 
longer than l~ , and becomes density-limited as soon as the MT length falls 
below i~ where the density profile is antenna-like. This implies that there is 
a discontinuous jump in the depolymerization speed right at L = t~ . 

(7) Finally, for all other values for St, the depolymerization speed of the MT re- 
mains rate-limited for lengths larger than a threshold length t&. At £4, which 
is smaller than i~ in this parameter regime, there is again a discontinuous 
jump to a density-limited depolymerization dynamics. 

If the depolymerization rate is larger or equal to the hopping rate of molecular 
motors, St > 1, then 5r > p^ a is always obeyed simply because pL a < 1. In this 
regime, all of the molecular details of the depolymerization kinetics are irrelevant. 
Neither cooperativity nor the actual value of the depolymerization rate matters 
in terms of the depolymerization speed; instead, only the bulk density regulates 
the speed. Note that this was the case for the data shown in Fig. [2j where we 
tentatively made the parameter choice St = 1. If the motors are faster than the 
depolymerization process, St < 1, we have to distinguish between the parameter 
regimes (a, (3, and 7, Fig. ^jfJ). Here the value of the depolymerization rate matters 
if the bulk density exceeds a certain threshold concentration, pL a > St, and the 
MTs are long enough. Finally, the depolymerization speed always becomes density- 
dependent and hence length-dependent if the MT length is short enough; the 
corresponding threshold length is £ reg = mm[£~ , id}- 

End-residence time strongly depends on cooperativity 

In contrast to the depolymerization speed, the mean end-residence time r res is 
strongly affected by the degree of cooperativity. Fig. [6] displays r res as obtained 
from our stochastic simulations for noncooperative and fully cooperative depoly- 
merization kinetics. Our simulations show that the end-residence time for the fully 
cooperative model is identical to the average lifetime of a terminal tubulin dimer 
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r res = T d '■= a /Kiepoi (Fig. [6]/l). Even for the noncooperative model, equals r d 
for large residence times and deviates from it only at small values. The relatively 
sharp transition to a constant lifetime of the terminal tubulin dimer occurs right 
at r r gg = r/pLa, i.e., the end-residence time equals the waiting time for a molecu- 
lar motor to arrive at the MT tip. For < r/pLa, the lifetime of the terminal 
tubulin dimer is identical to the arrival time (Fig. |6jv4 and B). Once the arrival 
time becomes shorter than the inverse depolymerization rate, the end-residence 
time levels off at t™ s = 1/5. These results show that the dependence of the end- 
residence time on density can be used to quantify the degree of cooperativity. 
This would require experiments with motor densities on the MT larger than those 
studied up to now (1T5| [23]) . 

The observation that the depolymerization speed is independent of the degree 
of cooperativity seems to be at odds with the experimental finding that the end- 
residence time, r res , of Kip3p depends on the total Kip3p concentration and is 
inversely proportional to the macroscopic depolymerization speed ( |23|) . Actually, 
however, there is no contradiction and the findings are readily explained within 
our theoretical model: For a noncooperative model, t^ s is simply given by the 
depolymerization rate, because after they arrive, the particles stay at the tip until 
they depolymerize the MT: 

C = ] ■ (14) 

For a fully cooperative model, r r f g S depends not only on 5, but also on the rate 
at which the second-to-last site becomes populated. Say the probability for the 
second-to-last site to be occupied is p + . Then, r r f g S is given by a sum of two 
contributions arising from the cases in which the second-to-last site is empty or 
occupied, respectively: 

T r f : s = (i-p + )(-j- + ^+ P+ ^. (is) 

If-the-second to last site is empty (which is the case with probability 1 — p + ) r res 
is the sum of arrival time r/pL a and depolymerization time 1/5. Otherwise, the 
end-residence time r res simply equals 1/5. 

As shown in the previous section, two distinct scenarios arise: For small bulk 
densities such that < 5t, the density profile at the plus-end exhibits a micro- 
scopic spike with p + = p\ ja ,/5r. For large densities, pL a > 5t, a macroscopic traffic 
jam emerges such that p + = 1. This result obtained for the motor density at the 
MT tip (Eq. [9]) may now be used to calculate r r f g S using Eq. 15 
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FIGURE 6: Motor end-residence times r res for cooperative and noncooperative 
depolymerization. {A) Mean end-residence time r rea plotted against the mean 
depolymerization time r^. Data were recorded for a range of depolymerization 
rates 5 = 0.02. . .2v. Noncooperative {shaded) and cooperative {black) dynamics 
are shown for different densities. {B) Mean end-residence time function 
of the Langmuir density pL a for various depolymerization rates (in units of v). 
For noncooperative depolymerization, r res is given by 1/5 {shaded lines). For the 
fully cooperative scenario {symbols), r res depends on whether the system is in the 
density-limited {St > pL a ) or in the rate-limited {St < pL a ) regime. While, for 
St > pLa, the end-residence time is given by r res = r/pLa {solid gray line), for 
St < pLa, h is density-independent and determined by the microscopic depolymer- 
ization rate r res = 1/5 (see also Eq. 16). 
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This agrees well with the results from stochastic simulations displayed in Fig. [6j A 
comparison with Eq. [6] shows that the end-residence time equals the typical depoly- 
merization time, i.e., the expected lifetime of a terminal tubulin dimer, r r f g S = r^. 
This is in agreement with experimental findings regarding the unbinding-rate of 
motors at the plus-end (l23j) and strongly supports the conclusion that depoly- 
merization of MTs by Kip3p is fully cooperative. Varga et al. (I23p measured the 
end- residence time of motors on double stabilized MTs, i.e., where depolymer- 
ization is switched off. They observed that the end-residence time is inversely 
correlated with the concentration of Kip3p, and fit their data with an exponential 
using a cut-off. This is in accordance with our results shown in Fig. [6]R However, 
since depolymerization has been switched off in the experiment, the rate 5, corre- 
sponding to the cutoff, now has to be interpreted as an unbinding-rate of motors 
at the plus-end. It would be highly interesting to design experiments where the 
depolymerization kinetics remains switched on, because this would allow one to 
measure the magnitude of the microscopic depolymerization rate S. 

DISCUSSION 

In this work, we have analyzed the effect of crowding and cooperativity on the 
depolymerization dynamics of microtubules. To that end, we constructed an 
individual-based model for the coupled dynamics of plus-end directed motor traffic 
and microtubule depolymerization kinetics. The model is based on well-established 
molecular properties of motors from the kinesin-8 family, i.e., the motors move on 
single protofilaments with high processivity at an average speed v, and exchange 
of motors between the bulk and the microtubule follows Langmuir kinetics. All 
parameters of the model, including the average walking speed, run length, and at- 
tachment rate, were directly extracted from available in vitro data (1231) . We have 
validated our model by reproducing the onset of length-dependent depolymeriza- 
tion as studied recently (1T51 1231) . Without using any additional fitting parameter, 
we found the same regimes of density profiles and ensuing depolymerization dy- 
namics as in the experiments, i.e., a linear antenna-profile with a length-dependent 
depolymerization speed and a flat profile with a constant depolymerization speed. 
Moreover, we identified a threshold density of motors above which a crowding- 
induced traffic jam emerges at the minus-end. The predicted shape and extent of 
these traffic jams should be amenable to experiments that raise the depolymerase 
concentration c or changing its rates of binding to and unbinding from the MT. 

The interplay between motor traffic and depolymerization kinetics at the mi- 
crotubule plus-end leads to strong correlations between the depolymerization dy- 
namics and density profiles of depolymerases. The plus-end acts as a bottleneck 
and crowding effects cause traffic jams. We find two qualitatively distinct regimes: 
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Motor densities below a critical threshold value, p* La = 5t, always show a local 
spike-like perturbation at the plus-end, the extent of which is the size of a het- 
erodimer. Above this threshold density, macroscopic traffic jams may emerge. 
These distinct density profiles at the plus-end affect the depolymerization speed 
and the end-residence time in qualitatively different ways. A quantitative analy- 
sis of the model using stochastic simulations as well as analytical calculations led 
to the following main results: The end-residence time of a depolymerase strongly 
depends on the degree of cooperativity. Whereas for noncooperative depolymeriza- 
tion kinetics the end-residence time is given by the microscopic depolymerization 
rate 5, it is density-dependent in the fully cooperative case: Increasing the Lang- 
muir density above the threshold value pj^ a = ^ r > the end-residence time changes 
from being inversely proportional to the density pLa to a constant value 5~ l . These 
results suggest an interesting way to determine the cooperativity of depolymeriza- 
tion kinetics and measure the value of the depolymerization rate 5. Although when 
the concentration c is increased, the end-residence time should be independent of 
concentration for noncooperative kinetics, it should strongly depend on concen- 
tration in the cooperative case. Experimental evidence points toward the latter 

In contrast, the depolymerization speed does not depend on the degree of coop- 
erativity of the depolymerization kinetics. Noncooperative and fully cooperative 
versions of the model give identical results. As a function of depolymerase concen- 
tration and the MT length, the depolymerization dynamics exhibits two qualita- 
tively distinct regimes: The depolymerization speed is either density-limited and 
determined by the bulk density of molecular motors, p-(x), or rate-limited and 
dictated by the value of the microscopic depolymerization rate, 5. Both regimes 
emerge due to crowding of molecular motors at the plus-end which acts as a bot- 
tleneck for molecular traffic. 

Density-limited regimes are correlated with microscopic traffic jams ("spikes") 
at the plus-end: The density profile self-organizes into a shape that cancels out all 
the effects of the depolymerization kinetics such that the depolymerization speed is 
solely determined by the bulk motor density, p_(x), and the average motor speed, 
v. Note that only in this regime length-dependent regulation is possible since the 
density changes over the MT length. As emphasized above, if the depolymerization 
rate 5 is larger than the hopping rate of the molecular motors, 5 > v, this remains 
the only regime of depolymerization dynamics. Then, the depolymerization speed 
is limited by the velocity of the plus-end directed motors, which is in accordance 
with recent experimental findings for Kip3p (123|) . In a parameter regime where 
motors depolymerize more slowly than they walk, 5 < v, there is a second rate- 
limited regime above the threshold density and for microtubules longer than 
some threshold length £ reg where Vd epo i = 0,5. In this regime the plus-end acts 
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as a strong bottleneck for molecular traffic. This causes a macroscopic traffic jam 
such that the motor density steeply rises to full occupation of all lattice sites at the 
plus-end of the microtubule. The cellular system sacrifices its capability to regulate 
the speed of depolymerization and only regains it once the MT length falls below 
£ reg , where the depolymerization speed again becomes density-regulated. From an 
evolutionary perspective one might speculate that the system has evolved towards 
8 = u, because this would allow regulation of the depolymerization dynamics over 
the broadest possible range. 

Beyond these observations, other predictions of our stochastic model can be 
put to test in experiments. By varying the motor concentration, two interest- 
ing observations could be made: First the phase diagram for the density profiles 
at the minus-end could be scrutinized experimentally. Second, the predictions 
on the density-profiles at the plus-end and their predicted strong correlations to 
the macroscopic depolymerization dynamics might be accessible to single-molecule 
studies. Manipulation of the molecular properties of the motor (e.g., the run 
length, attachment rate (H2|) . average speed and depolymerization rate) would 
change the intrinsic biochemical rates of the system and potentially lead to new 
parameter regimes. In addition, our results regarding length- and concentration 
dependence of the depolymerization process might be relevant in vivo, e.g., for 
mitotic chromosome alignment ( JT8j) . In our theoretical studies we explored the 
full parameter range, and therefore clear predictions are available for comparison. 

We believe that in a more general context, our theoretical work provides new 
conceptual insights into the role of collective and cooperative effects in micro- 
tubule assembly and disassembly dynamics. Future research could focus on the 
antagonism between polymerases and depolymerases (fT2 | H3 | Hill . spontaneous MT 
dynamics mediated by GTP-hydrolysis, the abundance of molecular motors in a 
cell, or more-detailed modeling of molecular motors (1261) . This may finally lead to 
a molecular understanding of the regulatory mechanisms of cellular processes in 
which MT dynamics plays a central role. 
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In this Supporting Material, details concerning the mathematical formulation 
and the stochastic simulations are given. In particular, the density profiles and 
the domain wall positions at the minus- and the plus-end are derived analytically. 
Further, some additional results are provided: (i) We show that the shapes of MT 
depolymerization curves (kymographs) are to a large extent independent of the 



choice of the depolymerization rate 5; see Fig. SI (ii) Analytical and numerical 



results from our theory are compared to experimental data on the relation between 



depolymerization speed and motor current ([!]); see Fig. S2 



Mathematical formulation 

In this article, we employ a lattice gas model. Its state is described by a set of 
occupation numbers n,i G {0,1} where i = 1,...,L denotes the lattice sites. In 
contrast to the notation in the main text, we here choose units of length and time 
such that the hopping rate from site to site and the lattice constant are both set 
to one. For an analytical description of the steady state density profiles of the 
molecular motors along the MT we consider the ensemble-averaged densities and 
currents: 

pi ■= (rii) , (SI) 
J i :=(n,(l-n i+1 )}. (S2) 

Note that the current Jj accounts for particle exclusion: a particle at site i moves 
to site i + 1 at rate v — 1 only if site % + 1 is unoccupied. The steady state results 
from a local balance between the transport current 

Jf := (rii_i(l - Hi)) - {rii{l - n i+l )) , (S3) 

the particle exchange with the bulk, 

Jl a := cw a (l - rii) - u d (m) , (S4) 
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and the depolymerization current, which sets the boundary condition at the plus- 
end. We now perform a mean-field approximation, where all spatial correlations 
are neglected, and a continuum limit keeping only the leading order terms ([3]). 
Then, the transport current simplifies to, 

J T (x) = (2p(x)-l)d x p(x), (S5) 

i.e. the transport current is proportional to the density gradient like a diffusion 
current in Fick's law but modified with a density-dependent prefactor which reflects 
site-exclusion between motors. The Langmuir current is given by 

J La (x) = cw a (l - p(x)) - uj d p(x) . (S6) 



Density profiles at the minus-end 

Within the above introduced framework the motor density profiles on the MT 
can be calculated analytically. In particular, the domain wall position, can be 
derived exactly as well as upon employing a linear approximation for the density 
profile close to the minus-end. For simplicity, we first consider the latter, especially 
because its results approximate the exact solution rather well over a broad range 
of parameters. 

Linear approximation 

In the immediate vicinity of the minus-end (x = 0) the density is small such that 
the full equation for the current balance, J T + J La = 0, 

(2p(x) - l)d x p{x) + cu a {l - p{x)) - u d p{x) = , (S7) 

reduces to d x p = cu a , which is solved by a linear (antenna) profile: 

P-(x) ~ cuj a x . (S8) 

At sufficiently large distances from the minus-end the density profile becomes flat. 
Therefore, J T vanishes and the system is dominated by the Langmuir kinetics, 



J = 0. Then, an asymptotic solution of Eq. S7 is given by the Langmuir density 



1 + K cu a + u d 

The full density profile is obtained by concatenating the antenna profile and the 
flat Langmuir profile such that the (local) current is continuous along the MT. 
There are two qualitatively distinct scenarios. For low bulk concentrations of 
molecular motors, c, the antenna profile matches the asymptotic Langmuir density 
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continuously resulting in a wedge-like profile; compare Fig. 2E in the main text. 
Approximately, the matching point, p_(d~) = pL a , is 

K 



V CM 



(SIO) 



cuj a (K+ 1) ' 

In contrast, above a certain threshold value for the bulk concentration, determined 
by K~ = 1, the two profiles can no longer be matched continuously and the 
density profile displays a localized discontinuity 02]), also termed a "domain wall" 
(DW). Its position is determined by a local current continuity condition fl2] [3]), 
p_(d~) = 1 — pLa, and can again be estimated using the linear antenna profile: 

1 



cu a {K + l)' 

Taken together, Eq. (5) from the main text is obtained, 

for K < 1 , 



(Sll) 



u *<- K , +r > ' (S12) 

1 for K > 1 . 



. CLU a (K+l) 

Exact solution and domain wall position 



To obtain the full solution Eq. [ST] has to be solved as already demonstrated in 
Ref. (]3]). Introducing a rescaled density at the minus-end cr_(s) = (2p(a;) — 1) — 



1 in Eq. S7 a transformed differential equation can be obtained 

(K + iy 



d x aJx) + d x In |cr(x)_| = u d - 



K 



which is mathematically equivalent to Eq. S7 and can be solved analytically 



a_(x) = iy_i(-K_(x)). (S13) 

Here W-\ is the second real branch of the Lambert iy-function (J3]) and Y_(x) 
reads (EI) 



Yjx) 



-IK 



K 



, (K + l) 2 2K 
exp <{ u S K _l x - j^-j \ . (S14) 



Herein the boundary condition p_(0) = corresponding to cr_(0) = —2K/ (K — 1) 
has already been accounted for. The local current condition for the domain wall, 
P-(d~) = 1 — PLa which corresponds to o"_(rf _ ) = —2 for the rescaled density, now 



enables us to calculate the DW position. Combining this condition with Eqs. |S13 
and IS 141 leads to 

d ( Ui , A ) = - i)2 . (S15) 
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Density profiles at the plus-end 



Analogously to the minus-end, we now evaluate the density profiles and ensuing 
DW at the plus-end. Because the tip steadily depolymerizes, the calculations have 
to be performed in a comoving frame which is introduced first. 

Comoving frame 

In the comoving frame, the above defined last lattice site L, i.e. the plus-end, 
is defined as the first site of the MT. This is equivalent to reverting the motor 
movement. Since molecular motors in the comoving frame move towards the first 
site of the lattice, the transport current changes sign J J = —jj' 



J^' ° := (n i+ i(l - in)) - (ni(l - n<_i)) , (S16) 

using the mean-field approximation as introduced above this leads to 

J T -> c °(x) = -(2p-l)d xP (x). (S17) 

The particle adsorption/desorption current jf a is unaffected. However, there is 
another contribution to the current balance in the comoving frame due to depoly- 
merization: Similar to the above definitions of the local currents, a local current 
which accounts for depolymerization in the comoving frame arises 

Jf° = 5(n i+1 - m) . (S18) 

Employing a mean- field approximation this expression simplifies to, 

j Co ( x ) = 5d xP {x) . (S19) 

This current term can be understood as follows. Due to the depolymerizing activity 
of a motor at the plus-end, in the comoving frame all motors on the MT simul- 
taneously approach the plus-end. In summary, by introducing a comoving frame 
the mean-field equation for the density at the MT plus-end p+(x) is obtained. In 
the steady state it reads 

(2p+ - 1 - S)d x p+ + cu a (l - p+) - cu dP+ = . (S20) 
Density profiles 

The above equation is solved in close analogy to Eq. |S7} In terms of a rescaled 
density 



2p + (x) 



K 
'K+l 



°+(x) = ' P - *+ , (S21) 



K 
'K+l 



[1 + 5) 



28 



the rescaled differential equation reads 

a' + (x) + d x In M*)| = k Ti K _ + { k1 1)5 ■ ( S22 ) 

The exact solutions to this equation are compared to stochastic simulations in the 
main text (Fig. 5^4). 

The solutions of this equation for 5 = are discussed in fl3]) and in parts above. 
However, in the case of depolymerization, i.e. for 5 > 0, two special solutions 
exist: Depending on the density of motors on the MT two classes of solutions for 
the density at the plus-end p+(x) can be distinguished. These are wedge-like or 
traffic jam density profiles; see Fig. 5 A in the main text. The boundary condition 
for these qualitatively distinct density profiles is p+{L) = 1. Defining 

Y(x) = M L )\^p{ K ^f K ^K 2 + l) {x ~ L) + a+{L) ] ' (S23) 
the two solutions for density profiles in the main text (black lines in Fig. 5^4) are 

(x) = {pLa + U 2 PLa - (1 + S))W (Y(x)) wedge-like 
\pLa + U 2 PLa-(i + S))W^(-Y(x)) traffic jam. 

Wo and W-i denote the first and the second real branch of the Lambert function. 
The reason for the form of these two solutions is the bottleneck (j5j) arising due to 
depolymerization (see main text). This bottleneck fixes the value of the tip density 
to its maximum p+(L) = 1 for pi a > 5. The transition from a traffic jam to a 
wedge-like density profile is thus not boundary-induced, i.e. due to a particular 
value of p+(L) < 1, but may be attributed to the depolymerizing activity of motors 
at the plus-end. As discussed in the main text, this transition is sharp and can be 
quantified in terms of p^ a and 5; see Fig. 5B. 

Domain wall position at the plus-end 

As already shown in the main text, microscopic jams can substantially influence 
the depolymerization dynamics. For large bulk concentrations, this perturbation 
no longer remains a local spike, but affects the profile on a macroscopic scale 
fl3J). Because the perturbation is macroscopic we can again use a hydrodynamic 
description, now with the boundary condition p+{L) = 1. Close to the plus-end 
Eq. S20 gives an approximately linear profile 



p+(x) w 1 - ^(L - x) . (S25) 

The slope increases with increasing depolymerization rate 5, c.f. Fig. 5A in the 
main text. In close analogy with the discussion for the minus-end there are two 
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scenarios for concatenating this linear profile with the Langmuir density. Here, for 
large enough Langmuir density and/or small enough depolymerization rates, we 
obtain a wedge-like profile with a matching point given by p + (d+) = pL a : 

compare the green curves in Fig. 5^4 in the main text. Upon increasing the depoly- 
merization rate or decreasing the Langmuir density a DW emerges whose position 
can be determined using current conservation, p + (d + ) = 1 — p^ a + 5: 



d + w L- 



K -6 



1 + K 



— ; (S27) 



compare the red curves in Fig. 5A The DW is most pronounced for K < 1. The 
height of the DW vanishes as K approaches the threshold value 

Kt = [±| (S28) 
from below. Equivalently, for a given K, the critical depolymerization rate reads 

4=^ = 2^-1='^. (S29) 
K + l cu a + u d 

Together with the condition for spikes, pLa < 8, this relation organizes the shapes 
of the density profiles into three classes: Microscopic jams at the tip, wedge profiles 
and DW profiles. 



Depolymerization dynamics of the antenna profile 

In the main text we have discussed how a spatially uniform density pL a affects 
depolymerization dynamics. Here we briefly show that our approach is also ap- 
plicable to linear antenna profiles, i.e. for MTs shorter than a certain threshold 
length, L < cf. Eq. (S12) and main text. Just as in the main text we equate 
the particle loss current due to depolymerization, 

</depol(X) = I 1 ~ P-( X )] P+( X ) S , (S30) 

and the particle flux towards the plus-end, 

J-{x) = [l-p-{x)]p-{x), (S31) 
and find 



p-0) 



for L < 



P+(x) = {, 5 fo^T' (S32) 



CU1„ 
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Since, according to Eq. (2) in the main text, the density at the plus-end deter- 
mines the depolymerization speed, Vd epo i — ^P+{ x )i the position-dependence of the 
tip density maps to a length-dependence of the polymerization speed. For MTs 
shorter than but longer than a certain depolymerization length l& = 5/cu a the 
depolymerization speed is length-independent 

Vdepoi = P+5 = 5 . (S33) 

Analogously to the result for constant bulk densities, this result shows that for MTs 
longer than the depolymerization length £& the dynamics of depolymerization of 
the antenna profile can not be distinguished from the dynamics as induced by a flat 
density profile. In contrast, at a MT length shorter than l& the depolymerization 
speed becomes length-dependent and follows the shape of the antenna density 
profile p-{x): 

V depol (L) = p + (L)5 = p^(L) . (S34) 

These results generalize the rate-limited and density-limited regimes discussed in 
the main text to non-uniform densities. Moreover, they show that once filaments 
become shorter than £_, i.e. the density profile is antenna-like, there is a second 
spike-induced length scale which is the relevant length scale for the onset of 
length-dependent depolymerization of MTs. 

Combining these results with the analogous conditions for the Langmuir plateau 
discussed in the main text, leads to the depolymerization regimes summarized in 
Fig. 5C and Table [TJ Simply put, the depolymerization dynamics changes from 
rate-limited to density-limited when the bulk density falls below the threshold den- 
sity 5: p~{L) < 5. For Langmuir densities below the threshold density, pL a < 8, 
the bulk density remains below the threshold density for the whole MT length 
such that the depolymerization dynamics is always density-limited and given by: 
Kiepoi — P-(L). This corresponds to regime (a) in Fig. 5C. For Langmuir densities 
above the threshold, pLa > 5, the depolymerization dynamics is rate-limited in 
the Langmuir plateau and given by Vdepoi — 8- in the antenna-like regime of the 
density profile, i.e. for L < £~, we have to distinguish between two cases: (i) K > 1 
(PLa > 0.5) where the bulk density profile exhibits a domain wall, and (ii) K < 1 
(PLa < 0.5) where the bulk density profile is wedge-like. In the latter case, the 
bulk-density profile changes slowly and hence remains above the threshold S for 
some time even below L = t~ . Only for MTs shorter than £4, given by P-(£d) = S, 
the dynamics changes from rate- to density-limited. This corresponds to regime 
(7) in Fig. 5C. In contrast, for K > 1 (pL a > 0.5), the bulk density p_(x) ex- 
hibits a discontinuous jump from the Langmuir density to 1 — pLa right at £~ . If 
this maximum value of the antenna-like profile is less than the threshold density, 
1 — PLa < 8, then the depolymerization dynamics discontinuously switches from 
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rate-limited to density-limited. This defines regime in Fig. 5C. Otherwise, if 
1 — PLa > S, we are back to regime (7). In summary, all regimes show a constant 
polymerization speed for long MTs in the Langmuir plateau. Depending on the 
relative magnitude of the Langmuir density and the depolymerization rate this 
regime is either density-limited and given by pLa or rate-limited and given by 5, cf. 
second column in Table [TJ In all scenarios the dynamics becomes length-dependent 
at some scale which is, however, different. While for regimes a and /3, it coincides 
with the beginning of the antenna-like density profile i~ , it is given by Id for regime 
(7) cf. third column in Table IT] 



Regime 


Condition 


Constant Vd epo i 


Critical MT length 


(a) 


PLa < 8 


Vdepoi = PLa 


l. 


(/?) 


PLa > S > 1 - p La 


Vdepoi = 5 


£_ 


(7) 


else 


Vdepol = 5 


id 



Table 1: Summary of the similarities and differences in the depolymerization 
regimes (a), (/3) and (7). While for MTs longer than a critical length (third col- 
umn) the depolymerization speed is constant (second column), it becomes length- 
dependent below this critical length and is then given by Vd ep oi — P-{L). 



Numerical implementation 

The stochastic dynamics of the individual-based model was simulated using a 
Gillespie algorithm ((7J) and employing the rates introduced above. Note that this 
method provides the mathematically exact stochastic dynamics. This is essential 
for the investigation of dynamic phenomena like length-dependent shortening. 

In Fig. 2A-D., our simulations of the motor traffic started from an initial con- 
dition where the MT lattice was empty and subsequently filled up with motors 
triggering the depolymerization dynamics. To visualize time- dependent MT length 
and motor densities in one kymograph we implemented 14 protofilaments and av- 
eraged the motor intensities and MT lengths, see Fig. 2A-C. In detail, the visual- 
ization of kymographs was achieved as described in the following. From stochastic 
simulation data of the MT, each second the occupation numbers of motors along 
the MT nf\ where j indexes the 14 protofilaments of the MT, were evaluated and 
converted to RGB color values: 

Ri = J_1 14 — — , Gi = 3 ~* 4 1 , Bi = 0. (S35) 

These values display the density of motors as green and the uncovered MT surface 
as red. Hence, if the MT is completely empty it is red, while at complete motor 
coverage it is green. 
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Steady state motor densities as shown in Figs. 2E and 5A were obtained by 
time-averaging over 10 4 independent realizations after an equilibration time of 
2000 time steps; note that for a constant lattice size time and ensemble averages 
yield identical results (J2J). In Fig. 5 A, the density profiles were recorded in the 
comoving frame of the MT plus-end, while density profiles in Fig. 2E were recorded 
in the rest frame of the MT minus-end. In Fig. 2E, the density profiles resulting 
from the minus end without any influences from the plus end are shown. This 
can be viewed as an infinitely long lattice. To simulate such a lattice, we chose 
the following boundary conditions. We neglected depolymerization as it arises 
at infinity. Further, we set the exiting rate for motors at the last site equal to 
1 — PLa- Then, the transport behavior at the tip is the same as on the lattice, if 
the Langmuir density is reached, p(x) = p^&- 

In the second part of the article we focus on the dependence of the depolymer- 
ization speed on the motor density. To this end, simulations were performed in 
a comoving frame where the density at the minus-end was fixed to the Langmuir 
density. This was achieved by extending the lattice one site to left with each de- 
polymerization step and filling the thereby created site with the probability p^ a . 
This procedure may also be interpreted as an infinite MT allowing to observe motor 
dynamics at the MT tip without perturbations arising from the length-dependent 
depolymerization regime. 

We measured the mean end-residence time of individual motors at the plus-end 
r res and the mean lifetime of the terminal tubulin dimer t&. Data of these were 
obtained by averaging over 10 4 time steps r, cf. Figs. 4 and 6. 



How kymographs become independent of the de- 
polymerization rate 

In Fig. [Si] we provide data that explicitly shows the parameter independence of 
MT depolymerization. This results has been generalized in the main text to all 
possible motor concentrations and depolymerization rates. 



Comparison with experiments: dependence of the 
depolymerization speed on the bulk flux and bulk 
density 

Experimentally it was found that the depolymerization speed is linearly correlated 
with the flux of molecular motors towards the plus-end (JT]). We have collected 



data from our simulations similar to experiments. Figure S2 shows a scatter plot 
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FIGURE SI: Kymographs as they become independent of the depolymerization 
rate 5. Time-space plots of depolymerizing MTs for different depolymerization 
rates are shown, ranging from 5 = 0.1 v (red) to 1.0 v (blue), the latter value cor- 
responds to the motor speed of v = 6.35 sites s . For slow depolymerization rate, 
S J$ 0.5, the depolymerization speed is related to the microscopic depolymeriza- 
tion rate 5, whereas for rates 5 ^ 0.5 the depolymerization speed is independent 
of the depolymerization rate but depends on the density of motors on the MT, as 
outlined in the main text. 



for the depolymerization speed as a function of the bulk flux of motors, J^, for 
two values of the microscopic depolymerization rate S. 

The noise in the ensemble of realizations has two sources. The bulk current 
fluctuates since the Langmuir kinetics responsible for the bulk density pLa is a 
stochastic process. The depolymerization speeds vary from realization to realiza- 
tion because the depolymerization kinetics is a Poisson-like process. Also shown in 



Fig. S2 are ensemble averages. These mean values, as predicted in the main text, 
show the following behavior. For a macroscopic depolymerization speed lower than 
the depolymerization rate, AL/AT < da, it is density-limited and identical to the 
bulk density: 

p^v for p La < 5t . (S36) 



AT 

Rewriting this relation in terms of the bulk current means that the data should 



fall on the parabola displayed as the solid curve in Fig. S2 For low densities 



PLa ~ 0.25, where crowding effects are weak, this implies AL/AT = p La t> « J La as 



observed experimentally (CE]); see Fig. S2 



As the bulk density is increased two things happen. First, crowding effects 
become important invalidating the linear relationship between bulk current and 



34 




FIGURE S2: Depolymerization speed AL/AT as a function of the bulk motor cur- 
rent on the MT, Ji a . Data from stochastic simulations for two different depoly- 
merization rates 5 = 0.5 v (o) and 0.8 v (<) are shown, each for a set of Langmuir 
densities (colors) corresponding to concentrations c = 0.16 . . . 13 nM. To illustrate 
the effect of statistics, small symbols show individual measurements as obtained 
from 10 3 fixed time measurements AT = 500r, while large symbols indicate their 
mean (AL/AT). Good agreement with experimental data (□) measured in the low 
density regime is found (JTJ). Here, the theoretical prediction given by Jl & (solid) 
is hardly discernible from a linear best fit to experiments (dashed). 



depolymerization speed. It would be interesting to test our prediction that the 
depolymerization speed is linear in the bulk density by using higher motor con- 
centrations or changed biochemical rates such that K becomes significantly larger 
than 1. 

Second, if pLa > the depolymerization speed becomes rate-limited: 

— ^ = a5 for p La > 5t . (S37) 

This puts an obvious upper bound on the depolymerization speed. It cannot 
become larger than the microscopic rate of depolymerization at the plus-end. If 
the depolymerization rate is larger than the hopping rate of the molecular motors, 
5 > u, the depolymerization speed is, for all possible values of the bulk density, 
strictly given by the bulk density. 
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